Identification and genetic characterization of a recently identified enterovirus C116 in China.

Enterovirus C116 (EV-C116) is a new member of the enterovirus C group which is closely associated with several infectious diseases. Although sporadic studies have detected EV-C116 in clinical samples worldwide, there is currently limited information available. In this study, two EV-C-positive fecal specimens were detected in apparently healthy children, which harbored low abundance, through meta-transcriptome sequencing. Based on the prototypes of several EV-Cs, two lineages were observed. Lineage 1 included many types that could not cause EV-like cytopathic effect in cell culture. Three genogroups of EV-C116 were divided in the maximum likelihood tree, and the two strains in this study (XZ2 and XZ113) formed two different lineages, suggesting that EV-C116 still diffuses worldwide. Obvious inter-type recombination events were observed in the XZ2 strain, with CVA22 identified as a minor donor. However, another strain (XZ113) underwent different recombination situations, highlighting the importance of recombination in the formation of EV-Cs biodiversity. The EV-C116 strains could propagate in rhabdomyosarcoma cell cultures at low titer; however, EV-like cytopathic effects were not observed. HEp-2, L20B, VERO, and 293T cell lines did not provide an appropriate environment for EV-C116 growth. These results challenge the traditional recognition of the uncultured nature of EV-C116 strains and explain the difficulty of clinical detection.

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts.

Coxsackievirus A16 (CVA16) is a major pathogen that causes hand, foot, and mouth disease (HFMD). The recombination form (RF) shifts and global transmission dynamics of CVA16 remain unknown. In this retrospective study, global sequences of CVA16 were retrieved from the GenBank database and analyzed using comprehensive phylogenetic inference, RF surveys, and population structure. A total of 1,663 sequences were collected, forming a 442-sequences dataset for VP1 coding region analysis and a 345-sequences dataset for RF identification. Based on the VP1 coding region used for serotyping, three genotypes (A, B, and D), two subgenotypes of genotype B (B1 and B2), and three clusters of subgenotype B1 (B1a, B1b, and B1c) were identified. Cluster B1b has dominated the global epidemics, B2 disappeared in 2000, and D is an emerging genotype dating back to August 2002. Globally, four oscillation phases of CVA16 evolution, with a peak in 2013, and three migration pathways were identified. Europe, China, and Japan have served as the seeds for the global transmission of CVA16. Based on the 3D coding region of the RFs, five clusters of RFs (RF-A to -E) were identified. The shift in RFs from RF-B and RF-C to RF-D was accompanied by a change in genotype from B2 to B1a and B1c and then to B1b. In conclusion, the evolution and population dynamics of CVA16, especially the coevolution of 3D and VP1 genes, revealed that genotype evolution and RF replacement were synergistic rather than stochastic.